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ABSTRACT 

The diffuse high-latitude Ha background is widely believed to be predomi- 
nantly the result of in-situ recombination of ionized hydrogen in the warm inter- 
stellar medium of the Galaxy. Instead, we show that both a substantial fraction 
of the diffuse high-latitude Ha intensity in regions dominated by Galactic cir- 
rus dust and much of the variance in the high-latitude Ha background are the 
result of scattering by interstellar dust of Ha photons originating elsewhere in 
the Galaxy. We provide an empirical relation, which relates the expected scat- 
tered Ha intensity to the IRAS 100 /xm diffuse background intensity, applicable 
to about 81% of the entire sky. The assumption commonly made in reductions 
of CMB observations, namely that the observed all-sky map of diffuse Ha light 
is a suitable template for Galactic free-free foreground emission is found to be in 
need of reexamination. 
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1. INTRODUCTION 

Efforts to observe the all-sky distribution of diffuse Ha emission (Dennison et al. 1998, 
VTSSQ Gaustad et al. 2001, SHASSA0 Reynolds et al. 2002, Haffner et al. 2003, WHAMfJ 
have resulted in detailed maps, which revealed the presence of diffuse, highly structured 
Ha radiation at high Galactic latitudes. Although questions about the inferred ionization 
sources for diffuse, ionized gas outside of H II regions in late-type galaxies have not been 



fully resolved (e.g. Ferguson et al. 1996; Hoopes & Walterbos 2000; Seon 2009), the high- 
latitude Ha emission in the Milky Way has generally been interpreted as resulting from 



in-situ emission in an extended, ionized warm interstellar medium (WIM) (Gaensler et al. 



2008 Haffner et al. 2009). This was supported by studies of the [S II]/Ha ratio in a limited 
number of selected regions ( |Reynolds|19"8^| [19881 [Haffner et al.|2003[|Madsen et al.|2006[ ) and 



with radiative transfer models by Wood & Reynolds (1999), who concluded that scattering 



of Ha photons originating elsewhere in the Galaxy by high-latitude dust was not a significant 



component of the Ha background. Consequently, an all-sky Ha map by Finkbeiner (2003) 



produced through the combination of the VTSS, SHASSA, and WHAM data, has been 
used extensively as a template for Galactic free-free emission in the analysis of Galactic 
foregrounds of the all-sky CMB observations by WMAP ( Reynolds fc Haffner||2005 Bennett 



et al. 2003; Dickinson et al. 2003 Gold et al. 2010), with the fundamental assumption that 



the extinction-corrected Ha intensity is directly proportional to the line-of-sight emission 
measure in the foreground Galactic interstellar medium. 

The question of which physical processes are responsible for high-latitude Ha radiation 
was reopened recently through the analysis of the Ha excess observed in the isolated, high- 
latitude interstellar cloud LDN 1780 (1 = 359°. 02; b = +36°. 77; distance ~ 110 pc; |Franco 



(1989)). LDN 1780 is readily identifiable as a distinct object on the map of Finkbeiner (2003) 
by its excess Ha surface brightness. From an analysis of the Ha morphology, in addition to 



CO, FIR, and optical extinction maps of LDN 1780, del Burgo & Cambresy (2006) concluded 



that the observed surface brightness distribution of LDN 1780 was inconsistent with photo- 
ionization as the cause for the observed Ha light. Photoionization by an external field of 
Lyman continuum photons is expected to result in ionization of a thin outer shell of an 
isolated interstellar cloud while the observations of the Ha light from LDN 1780 clearly 
reveal the dense inner core of the cloud as the brightest feature (see Fig. 1 (bottom)), 
showing that most of the observed Ha photons come from the densest interior structures of 
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the cloud. To account for this phenomenon, del Burgo & Cambresy (2006) proposed that 



cosmic rays caused partial ionization of the hydrogen throughout LDN 1780, leading to the 
observed Ha emission. To produce the observed intensity, they postulated a locally enhanced 
cosmic ray flux about one order of magnitude larger than commonly accepted for the solar 
neighborhood. 



Mattila et al. (2007), analyzing the same LDN 1780 data from the Finkbeiner (2003) 



map as del Burgo & Cambresy (2006), concluded instead that the Ha excess intensity in 
LDN 1780 and a number of other similar high-latitude clouds is more likely the result of 
scattering by the dust in these clouds of Ha photons produced elsewhere in the Galaxy. Their 
interpretation relies upon the facts that (1) the optical depth of LDN 1780 at the wavelength 
of Ha is sufficiently small so that external Ha photons from the Galactic interstellar radiation 
field (ISRF) can readily reach the inner structures of LDN 1780 and (2) the known albedo of 
interstellar dust is sufficiently high to produce the observed intensities by scattering, given 
the known mean intensity of the Galactic radiation field in the Ha line. Still more recently, 



Lehtinen et al. (2010) analyzed the excess Ha surface brightness of the large high-latitude 



translucent cloud G294-24 and concluded that also there the excess Ha intensity can be 
explained solely by galactic Ha radiation scattered off dust in the cloud, thus supporting the 



conclusions of Mattila et al. (2007) 



The possible presence of scattered Ha radiation at high latitudes was first considered by 



Jura (1979), but subsequent Monte Carlo radiative transfer simulations by Wood & Reynolds 



(1999) suggested that the scattered Ha component should be relatively small. They conclude 



that the scattered light from midplane H II regions is typically around 10% of the total Ha 
intensity at high latitudes, with a similar fraction (see their Fig. 3) being contributed 
by scattering of Ha photons emitted by the WIM. Thus, their models predict an average 
scattered Ha contribution at high latitudes of about 20% of the total Ha intensity. Obvious 
shortfalls of the Wood & Reynolds model are their assumption of a smooth, unstructured 



ISM, and more seriously, the use of the Henyey & Greenstein ( 1941 ) scattering phase function 



The Henyey-Greenstein phase function is known to under-predict scattered light intensities 



e.g. Witt 1977; Draine 2003), when large scattering angles are involved as they are in the 



scattering of midplane Ha photons by high-latitude clouds. 

Frequently cited support for the conclusion that scattering makes only a minor con- 
tribution to the diffuse Ha background comes from WHAM measurements of the [S II] /Ha 
ratio cited earlier in this section, which demonstrated that this ratio is typically several times 
larger in the diffuse ionized interstellar medium than in bright H II regions. If scattering 
of Ha photons originating in H II regions were to dominate the Ha intensity in the diffuse 
ISM, the [S II] /Ha ratio would be expected to be similar to that in H II regions. We will 
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show in this paper that the selection criteria in use for the [S II] /Ha ratio measurements 
conducted with WHAM have led to a highly biased sample of regions in which scattering 
indeed is a minor contributor. The [S II] /Ha argument, while valid in principle, has been 
inappropriately extended to the entire sky without sufficient observational justification. 



Given the ubiquity of high-latitude dust structures in the Galactic cirrus (Sandage 



1976 Low et al. 1984), the conclusions of Mattila et al. (2007) and Lehtinen et al. (2010) 



have profound consequences. If a substantial fraction of the high-latitude Ha background is 
due to scattering rather than in-situ recombination, be that in balance to photo-ionization 



(Reynolds et al. 1995) or cosmic ray ionization as proposed by del Burgo & Cambresy 



(2006), the template used for Galactic free-free emission in the WMAP analyses (Bennett 



et al. 2003 Gold et al. 2010) employed systematically intensities that were too high, as well 



as spatial brightness variations not actually present in the distribution of the in-situ emission 
component alone. 

In this paper we present new observations of LDN 1780 in Ha and in adjacent continuum 



bands with superior spatial resolution and sensitivity compared to the Finkbeiner (2003) 
data. We use these new data to independently examine the case for dust scattering in 
LDN 1780. Scattering by dust would clearly be the phenomenon with greater consequences, 
because it does not require special conditions such as the locally, ten-fold enhanced cosmic 



ray flux postulated by del Burgo & Cambresy (2006). With the case for scattering made by 



Mattila et al. (2007) and by Lehtinen et al. (2010) confirmed, we then extend our conclusions 



to the entire high-latitude sky with consequences for future CMB analyses. We describe our 
new observations in Sect. 2, and discuss the calibrations and reductions in Sect. 3. In Sect. 4 
we analyze the excess Ha intensity found in LDN 1780. Sect. 5 contains the discussion of the 
consequences of our results for the high-latitude Ha background in general. Our conclusions 
are summarized in Sect. 6. 



2. OBSERVATIONS 

We used two identical, remotely operated Takahashi Epsilon 180ED astrographs (f/2.8) 
to obtain new imaging data for LDN 1780. The telescopes were located at New Mexico 
Skies Observatory near Mayhill, NM, at an altitude of 7300 ft (2225 m). Each telescope was 
equipped with a Santa Barbara Imaging Group STL-6303E CCD camera, covering a field 
of 2° x 3°. The angular resolution determined by the 9 jum pixels was 3.70" px _1 . Each 
camera was equipped with an 8-position filter wheel, jointly capable of housing 15 BATC 



intermediate-band filters (Yan et al. 2000) plus a narrow-band Astrodon Ha filter. The 



filter characteristics, total exposure times, and observing dates are summarized in Table 1. 
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The data were part of a larger program directed at studying the optical spectral energy 
distributions of high-latitude clouds in conjunction with Spitzer infrared band observations 
with both the MIPS and IRAC imaging systems (Witt et al. 2010, in preparation). 

The observations were carried out under photometric conditions. Both astrograph/ camera 
systems were flux-calibrated through observations of the absolute flux standard star SAO 
90153. For LDN 1780, individual sub-exposures of 900 s each were processed through bias 
subtraction, flat-fielding and cosmic ray removal. Atmospheric extinction corrections were 
applied to each 900 s sub-exposure, before co-adding them. This was particularly important 
in the case of LDN 1780, because its negative declination required data acquisition over a 



fairly large range of air masses. We note here that the SHASSA data (Gaustad et al. 2001) 



which form the basis for the Finkbeiner (2003) Ha map for LDN 1780 were not corrected 
for atmospheric extinction. 

The narrow-band Ha filter has a nearly rectangular transmission profile with 6 nm 
FWHM. In addition to the Ha line, this filter also transmits scattered-light continuum 
around Ha and possibly light from the [N II] line at 654.81 nm and to a lesser extent light 
from the [N II] line at 658.34 line, if there is significant radiation present in these lines in 



LDN 1780. As discussed by Madsen et al. (2006), these lines are fairly weak relative to Ha in 



H II regions but are more significant in the radiation from the WIM. Observations of stellar 
continuum sources with spectral energy distributions similar to that of LDN 1780 with the 
BATC #9 filter and the Ha filter were used to determine the continuum contribution passing 
through the finite width (6 nm) of the Ha filter. 

On 13 July and 14 July, 2010, one of us (GJM) used the Wisconsin H-Alpha Mapper 
(WHAM), currently located at CTIO, to measure the emission line ratios [N II]/Ha for the 
stronger [N II] line at 658.34 nm and [S II]/Ha for the [S II] line at 673.1 nm on LDN 1780 
relative to the nearby sky. The 1-degree WHAM bearrj^] encloses the full extent of LDN 
1780 as well as as substantial fraction of the adjacent sky. The excess Ha intensity of LDN 
1780, averaged over the 1-degree beam relative to three sky positions at several degrees offset 
(same latitude, 1°.52 to lower 1, 1°.30 to higher 1; same longitude, 3°.0 to higher b) was found 
to be 0.90 R; the Ha line profile was substantially narrower (18 km/s) than commonly seen 
in the high-latitude WIM (25 - 30 km/s). The line ratios were measured as [N II]/Ha = 0.19 
and [S II] /Ha = 0.16, rather close to corresponding values in H II regions and much smaller 



than typically found in the diffuse WIM (Madsen et al. 2006). 



While these recent WHAM measurements suggest that significant contamination of our 
Ha intensities on LDN 1780 by contributions from the [N II] lines is not likely, we will 
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carry out an additional examination of this question in Sect. 3 by comparing our measured 
intensities with those derived from the SHASSA data, which have been calibrated with pure 
WHAM Ha intensities and should therefore be free of [N II] contamination. In the absence 
of contamination, we expect the two sets of Ha intensities to agree within the measurement 
uncertainties. 



3. ADDITIONAL DATA AND REDUCTIONS 

3.1. Surface Brightness Morphology 

We conducted differential surface brightness measurements for 378 individual circular 
source regions on LDN 1780 (radius 37") and the adjacent sky, represented by the average 
of 40 circular sky fields, as displayed in Fig. l(top). The positions of the source regions were 
selected by eye to avoid contributions from stars as faint as V ~ 21 mag, yet were aimed to 
provide maximum surface coverage of the visible extent of LDN 1780. 

The source regions were divided into five color-coded groups to distinguish regions by the 
dust column density through the cloud and the ISRF impinging upon the dust in them (Fig. 
l(top)). Fig. 1 (center) displays the image of LDN 1780 in the narrow-band Ha filter. The 
morphology of LDN 1780 observed with the Ha narrow-band filter appears to be identical 
to that in adjacent continuum filters. Given the width (6 nm) of the Ha narrow-band filter, 
both Ha line photons as well as continuum photons are contributing to this image. Fig. 
1 (bottom) displays the residual Ha image of LDN 1780, after the data in Fig. 1 (center) 
have been sky-subtracted and corrected for the continuum contribution in the narrow-band 
Ha filter. The latter was determined to be 12.38% of the continuum intensity measured 
through the broader BATC #9 filter, normalized to the same total exposure time, sky- 
background and extinction corrected. This image shows the excess Ha of LDN 1780 only. In 
order to overcome the large per pixel noise but still leave the stars as recognizable discrete 
sources, the image has been smoothed by a 7-pixel Gaussian. The most important aspect 



of the resulting image is that the high-opacity inner core region (Ay ~ 4 mag; Ridderstad 



et al. (2006)) of LDN 1780 appears as a region of greatest Ha brightness. The implications 



of this fact will be discussed below in Sect. 4. 



3.2. External Data 



Other data sets employed in this study were the all-sky Ha map of Finkbeiner (2003 ) and 
the all-sky maps of the IRAS 100 /mi diffuse background and the line-of-sight dust extinction 
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maps of Schlegel et al. ( 1998 ). Finally, we used numerical data for the optical depths through 
LDN 1780 that were derived from the intensity at 160 /im (Spitzer, MIPS) and the dust 
temperatures determined by Witt et al. (2010, in preparation). The dust temperatures 
were determined from the intensity ratios at 160 /im (MIPS) and 100 /xhi (IRAS) and thus 
refer to the equilibrium temperature of the large-grain population responsible for the far- 
infrared emission as well as the optical scattering in LDN 1780. The temperatures of these 
grains were found to lie in the range 14.5 K to 16.8 K, typical for dust in interstellar clouds 
exposed to the ISRF near the sun (Bot et al. 2009) and in good agreement with the dust 
temperatures for LDN 1780 independently determined by Ridderstad et al. (2006). We note 
that these dust temperatures are well below the minimum temperature (17.35 K) indicated 
by the Schlegel et al. (1998) dust temperature map for the center of LDN 1780. The reason 
for this discrepancy most likely arises from the low spatial resolution of the COBE/DIRBE 
100 /im and 240 /im maps, which are the basis for the Schlegel et al. temperature map. 
With our temperatures based on the much higher IRAS 100 /tm resolution, we determined 
the optical depth at 160 /im for each sampling position on LDN 1780 from the relation 
r i60 — Iwo/ BiQ (T), where B 1&Q {T) is the Planck function for the dust temperature at the 
specific sampling field and Ji 60 is the corresponding measured MIPS intensity, both at the 
wavelength of 160 /im. 



3.3. Possibility of Emission Line Contamination 



We examined our measurements of the continuum-corrected excess Ha intensities with 
regard to the reliability of our absolute calibration and the possible presence of [N II] line 
contamination by comparing them with the corresponding measurements of the excess in- 
tensities in the Finkbeiner Ha map. The two data sets differ in several important aspects. 
(1) Resolution: Our data (NMS, for New Mexico Skies) have a scale of 3.70" px _1 while the 
SHASSA survey, which underlies this region of the Finkbeiner Ha map, was conducted with 



an angular scale of 47.64" px 1 , with the final resolution of the Finkbeiner (2003) map re- 



duced to a resolution of 6'; (2) Star subtraction: Our data are free of starlight contamination 
to the limit of V ~ 21 mag by design while the SHASSA survey relied upon star subtraction 
based on stellar images in adjacent continuum images, which was not perfect in all cases 
(see Gaustad et al. (2001) and Finkbeiner (2003) for details); (3) Atmospheric extinction 



correction: The NMS data were corrected for atmospheric extinction at the sub-exposure 
(900 s) level while the SHASSA data where not corrected for atmospheric extinction; and 
(4) Absolute calibration: The NMS data were calibrated against the absolute flux standard 
star SAO 90153, whereas the SHASSA data were calibrated against sources with known Ha 
fluxes, with an absolute zero point based on WHAM Ha measurements, which are pure, 
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uncontaminated Ha intensities. The LDN 1780 data obtained from the Finkbeiner map 
should, therefore, be free of contamination from [N II] lines, while the NMS data may not 
be. 

The comparison of the intensities measured in LDN 1780 is shown in Fig. 2. In order 
to facilitate this comparison, the resolution of the Finkbeiner map was artificially increased 
to match that of the NMS data. The ideal 1:1 correlation is indicated by the solid diagonal 
line; the least-squares linear fit (dashed line) through the origin suggests that the Finkbeiner 
intensities, on average, are about 11% higher than the corresponding NMS intensities. This 
difference is small compared to the typical errors in each measurement represented by the 
error bars in Fig. 2. This independent test is consistent with our recent measurement of the 
[N II] /Ha line intensity ratio for LDN 1780 (see Sect. 2), which showed the [N II] lines to 
be weak. We, therefore, conclude that contamination of the NMS data by light from [N II] 
lines passing through our Ha filter is insignificant. 



ANALYSIS OF THE Ha EXCESS IN LDN 1780 



As seen in Fig. 1 (bottom), LDN 1780 stands out brightly against the surrounding sky at 
Ha. We further observe that both sets of measurements shown in Fig. 2 consistently identify 
the densest parts of LDN 1780 with Ay ~ 4 mag, represented by the red symbols (see Fig. 
l(top)), as the region of maximum Ha intensities. Such brightness contrast would not be 
expected, if the Ha radiation from LDN 1780 were the result of photo-ionization. If ionizing 
photons could penetrate via low-density paths and voids through a turbulent interstellar 



medium to the altitude of LDN 1780, as suggested recently by Wood et al. (2010), the 



resulting ionization would be limited to a thin outer shell of LDN 1780, given its maximum 
visual extinction of Ay ~ 4 mag and its known nature as a molecular cloud with a dense 



core and a more diffuse outer envelope (Toth et al. 1995). In such a morphology, the outer 



optically thin regions would appear brightest in Ha line emission, while the optically thick 
core regions would appear only about half as bright, as the extinction in the core would 
attenuate the Ha radiation from the far side. This is contrary to the observations shown in 
Fig. 1 (bottom). The Ha morphology of LDN 1780 is indistinguishable from its appearance 
in the adjacent continuum, which is well understood as originating through scattering of 



the ISRF by the dust grains in LDN 1780 (Mattila 1979). These facts strongly imply that 



photoionization by external Lyman continuum photons is not the cause of the observed Ha 
light in LDN 1780, while dust scattering of Ha photon from elsewhere in the Galaxy could 
provide a natural explanation. 



The same reasoning caused 



del Burgo & Cambresy 



(2006) to look for an ionization 
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source other than photoionization for LDN 1780, one that would be capable to penetrate 
deeply into the interior of the cloud. They suggested ionization by cosmic rays, which would 
be supported by a linear relation between the resulting Ha intensity and the optical depth 
through the cloud. In our view, they were misled into believing that there exists such a 
linear relation between the Ha intensity and the visual extinction through LDN 1780 by 
applying what we consider an inappropriate extinction correction (their equation (7)). We 
will show in this section that the actual observed functional dependence of the Ha surface 
brightness on the optical depths of different cloud regions is a natural consequence of the 
radiative transfer of an external ISRF through the interior of LDN 1780. 

The actual observed scattering response of LDN 1780 to illumination by an external 
ISRF in terms of the intensity / as a function of the line-of-sight optical depth r is the 
solution to the equation of transfer. In the range of optical depths from the optically thin 
limit to moderately optically thick regimes as represented by LDN 1780, this solution has 
the general functional form: 



J(r) = A[l - exp(-r)] (1) 

The constant A is a product of three quantities: (1) the average intensity J of the 
incident ISRF; (2) the dust albedo a; and (3) the phase function contribution factor P, 
which measures the coupling between the anisotropic external ISRF to the line of sight to 



the cloud through the forward-directed scattering phase function of the dust grains (Witt 



1985). Rather than trying to find the solution to the actual transfer problem for LDN 1780 
through modeling with all its inherent approximations and uncertainties, we will employ an 
empirical approach by using the well-observed scattered light intensities in our two continuum 
filters immediately adjacent to Ha, BATC #8 (607 nm) and BATC #10 (705 nm), to provide 
the functional parameters for the solution to the equation of transfer. 

In Fig. 3a we plot the averaged observed intensities I cont in these two filters against the 
optical depth r 160 //m (see Sect. 3) for all observed positions on LDN 1780. The least-squares 
curve fit shown corresponds to the function: 



Icont = A[l - exp(-BT W0 )], (2) 

where the constant A is the appropriate value for the continuum adjacent to Ha and the 
constant B is the ratio of the optical depths around Ha and at 160 /im. We find these values: 
A = 31.86 x 10~ 9 erg cm -2 s -1 A~ x sr _1 and B = 1157.92. This lets us conclude that the 
maximum optical depth of LDN 1780 near Ha is = 0.0026 x 1157.92 = 3.01, consistent 
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with the maximum extinction through LDN 1780 of Ay = 4.2 mag found by |Ridderstad 



et al.| Q2006fl and the ratio of r Ha /A v = 0.714 QWeingartner fc Draine||2001[ ) for Milky Way 
dust with Ry = 3.1. The average continuum ISRF intensity in the vicinity of the sun near 
Ho; is estimated to be J(656 nm) = 1.5 x 10 -7 erg cm -2 s _1 A -1 sr _1 (Porter & Strong 2005) 
and the dust albedo is modeled as a = 0.56 (Draine 2003). This allows us to estimate the 
phase function contribution factor P = 0.38. A value of P less than unity is expected for 
forward-scattering grains and a scattering geometry where large-angle scattering prevails, 
as is the case of high-latitude clouds illuminated preferentially by sources near the Galactic 



(Witt 


1985 


Stark 


1995 



clearly non-isotropic. The southern regions facing the Milky Way plane, represented by blue, 
cyan and green symbols are illuminated by a stronger radiation field, whereas the northern 
regions represented by the yellow symbols are illuminated by intensities that are 20% to 30% 
lower compared to the southern regions. 

A relation identical in functional form to Eq. (2) should also fit the measurements of 
the Ha excess surface brightness, if the Ha light in LDN 1780 is a result of scattering by 
the same grains with the same scattering geometry as the ones that apply to the adjacent 
continuum. To arrive at the appropriate value for the constant A, the first factor representing 
the average continuum radiation density needs to be replaced by the average Ha intensity. 



This is found by averaging the Finkbeiner Ha map over the entire sky. Mattila et al. (2007) 



and Lehtinen et al. (2010) report a value for the all-sky Ha average intensity of 8 R. The 
other two factors contributing to the value of the constant A remain approximately the same 
as found earlier for the continuum, given that the illumination geometries for the continuum 
and Ha are very similar. Thus, we estimate A(Ha) = 1.70 R, if we assume that the ratio 
of the Ha line intensity to continuum intensity as seen by LDN 1780 is the same as found 
near the sun. The constant B should be identical to that found from Fig. 3a, because, by 
design, B refers to an effective wavelength of 656 nm, the same as Ha. 

We plot our Ha surface brightness measurements as a function of the optical depth 
r 160 in Fig. 3b. The solid curve is the expected scattering solution for A = 1.70 R. It fits 
the Ha intensities for the northern portion of LDN 1780 rather well but fails to predict the 
higher intensities in the southern and central portions of the cloud. A more satisfactory 
solution for these portions of LDN 1780 with A = 2.35 is represented by the dashed curve, 
which implies that the radiation field as seen by the southern side of LDN 1780 has a 38% 
higher ratio of Ha to continuum intensities than the ISRF near the sun. An examination 
of the environment of LDN 1780 on the Finkbeiner Ha map reveals two giant H II regions 
at Galactic latitudes b ~ +22° just to the south of LDN 1780, with maximum Ha surface 
brightnesses near 100 R. These two H II regions are powered by the OB stars £ Oph and 



5 Sco, with approximate Hipparcos-derived distances (Perryman et al. 1997) of 140 pc and 
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123 pc, respectively. These two H II regions are therefore separated from LDN 1780 by only 
about 50 pc to 40 pc to the south, illuminating LDN 1780 under albeit rather inefficient 
scattering angles near 70°. The presence of these H II regions could well explain the greater 
relative strength of Ha in the southern portions of LDN 1780. Considering the closeness 
of the fits in Fig. 3b, we conclude, in agreement with Mattila et al. (2007), that both the 
absolute Ha surface brightness in LDN 1780 and in particular its dependence on the optical 
depth through this cloud are fully consistent with dust scattering by typical interstellar 
grains in this cloud. These fits do not require a ten-fold locally-enhanced cosmic ray flux, 
nor are they consistent with the expected surface brightness distribution produced through 
photoionization by external Lyman continuum photons. 



5. CONSEQUENCES FOR THE HIGH-LATITUDE Ha BACKGROUND 



The high-latitude Galactic sky is dominated by extended dust structures, which have 
been observed at both optical ( Sandage]|1976 Witt et aL]|2008 ) and far-infrared wavelengths 
( Low et aL]|1984 ). The dust morphologies in the two spectral ranges mirror each other closely 
( de Vries fc Le Poole|1985 Stark|199~5 ), supporting the conclusion that the optical continuum 
intensity is due to scattering of the Galactic ISRF by wavelength-sized dust grains, while the 
far-infrared emission is due to the equilibrium thermal emission from the same dust grains, 



re-radiating that portion of the ISRF that is absorbed by them (Guhathakurta & Tyson 



1989 Guhathakurta & Cutri 1994). Both the absolute and relative intensities of these two 



radiation phenomena are consistent with the intensity of the Galactic ISRF in the solar 



vicinity and the known scattering parameters (Draine 2003) of interstellar dust. 



At high latitudes, the intensities of the optical continuum and far-infrared backgrounds 
are linearly related, because the respective intensities result from scattering and absorption 
of the same interstellar radiation field by the same population of interstellar grains, and 
because the high-latitude sky is optically thin at both optical and far-infrared wavelengths. 
However, only at blue wavelengths can the optical background be fully attributed to scat- 



tering (Mathis 1973), while the optical background at red wavelengths appears to exhibit a 



significant contribution from dust luminescence or extended red emission (ERE), as shown 



by Guhathakurta & Tyson (1989) and Gordon et al. (1998). Therefore, in the subsequent 



analysis we rely only on the observed relation between the intensities at blue optical and at 
far-infrared wavelengths. 



For optically thin cirrus clouds, Witt et al. (2008) found a linear relation between the 



optical scattered light intensities near 457 nm and the 100 /im thermal emission intensities 



recorded by the IRAS survey (Schlegel et al. 1998): 



7(457 nm) [erg cm" V^^sr" 1 ] = 3.243 x 1(T 9 7(100 fim) [MJy sr -1 ] (3) 
The optically thin high-latitude dust clouds exhibit a scattered light spectrum with an 



intensity ratio (Witt et al. 2008) 



7(656 nm) 
7(457 nm) 



0.68 



(4) 



From the results shown in Fig. 3a,b and the associated curve fits, we adopt a line-to- 
continuum ratio 



J(77a)/J(656 nm) = 6.0 x 10 7 [R erg" 1 cm 2 s 1 A 1 sr 1 ]. (5) 

This corresponds to a fit to the Ha intensities in Fig. 3b with a constant A = 1.92, 
which is intermediate to the two fits shown with A = 1.7 and A = 2.35 but closer to the 
lower one. 

Combining the last three relations, we expect a scattered Ha intensity in the high- 
latitude, optically thin dust corresponding to the measured far-infrared intensities 7(100 
\im) [MJy sr -1 ]: 



7 (77a) [R] = (0.129 ± 0.015) x 7(100 /xm) [MJy sr" 



(6) 



This relationship assumes that all high-latitude dust structures see an ISRF intensity 
approximately the same as that observed in the solar vicinity, with the same Ha line to 
continuum ratio. We justify this assumption with the fact that virtually all high-latitude 
cirrus dust clouds owe their large angular extent and their high latitudes to their proximity 



to the sun (Weiland et al. 1986 Franco 1989 Szomoru &; Guhathakurta 1999). The cloud 



distances, typically 100 pc to 200 pc, are small in comparison with the distances of most of 
the principal Ha sources at low Galactic latitudes. 



By applying Eq. (6) to the 100 /im IRAS map of Schlegel et al. (1998), we now have a 
way for estimating the contribution to the high-latitude Ha background due to Galactic dust 
scattering. We estimate the uncertainty in this prediction to be about 12%, a combination 
of the estimated uncertainties in the relations shown in Eqs. (3), (4), and (5) of 3%, 5%, 
and 10%, respectively. 



We applied Eq. (6) to the entire high-latitude sky, where the linear, single-scattering 
approximation Eq.(6) is expected to be valid (r(77a) < 1). Fig. 4a shows the resulting 
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histogram of the numbers of 6' pixels with measured total Ha intensities (black), predicted 
scattered Ha intensities (green) and residual Ha intensities (red) attributed to in-situ re- 
combination in the WIM. The most likely value of scattered Ha intensities (0.1 R) is about 
19% of the most frequently encountered measured total Ha intensity (0.52 R). Thus, when 
averaged over the entire high-latitude sky, our prediction is in fair agreement with the expec- 



tations from radiative transfer models of Wood & Reynolds (1999) and the recent analysis of 



the free-free emission and associated Ha emission from the WIM by Dong & Draine (2010). 



We next analyzed a northern Galactic sky region associated with extended cirrus struc- 
tures, limited by Galactic longitudes (340° to 60°) and Galactic latitudes (+35° to +75°). 
Fig. 4b shows the corresponding histograms of the frequency of the observed Ha intensities 



from the map of Finkbeiner (2003), of the predicted Ha scattered light intensities, and the 



residual Ha intensities that can be attributed to in-situ recombination, all without extinc- 
tion corrections for this region. Here, the total observed Ha intensity displays a distinct 
double-peaked distribution, while the predicted scattered Ha intensities peak at a higher 
value (0.16 R) than for the high-latitude sky overall (0.10 R). We find that the removal 
of the scattered Ha contribution leads to a narrowing of the frequency distribution of the 
residual Ha intensities attributable to in-situ emission, implying smaller amplitude in the 
spatial structure of the in-situ emission component than previously thought ( Hill et al.|2007 ), 
and an overall reduction in average intensity by about a factor of two. This implies that 
the free-free emission template for such high-latitude regions used in the WMAP reductions 



(Bennett et al. 2003 Gold et al. 2010) has been overestimated in both strength and spatial 



structure. We suggest that the standard assumption in the reduction of CMB observations, 



namely that the observed Ha all-sky map of Finkbeiner (2003) is a suitable template for 



free-free emission in the Galactic foreground, must be reexamined, especially as it applies to 
high-latitude parts of the sky. 

In contrast to the effect of scattering, the correction for foreground extinction, using the 



standard approach of Bennett et al. (2003), of the observed Ha intensities or the Ha residual 



intensities after removal of the scattered component introduces relatively small changes, as 
illustrated in Fig. 4c. 

The observed Ha high-latitude background before correction for scattered Ha is much 
more structured than the component we identify as attributable to in-situ recombination. 
The FWHM values of the two respective intensity frequency distributions shown in Fig. 4c 
are different by about a factor of two. Fig. 5 shows a map of all high-latitude parts of 
the sky with 1 degree resolution, where the linear, single-scattering approximation Eq. (6) is 
expected to apply (r(Ha) < 1). We display the ratio of scattered Ha to total observed Ha 
intensities. We find many extended regions where the fraction of the scattered Ha intensity 
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is of order 0.5 or higher, while it is much lower in other regions. Specifically, at a resolution 
of 1 degree, we find for the unmasked areas, representing 81.4% of the entire sky, that 51.7% 
exhibits a ratio of scattered-to-total Ha of < 0.25, 38.7% of the region falls into the range 
0.25 - 0.50, 8.3% reaches values of 0.5 - 0.75, with only 1.3% exceeding the value of 0.75. 
The last category increases in coverage about fourfold when the pixel size is changed from 
1 degree to 6', with a corresponding decrease in the second interval (0.25 to 0.50). This 
emphasizes the fact that the scattered Ha contribution adds particularly to the small-scale 
structure of the Ha background. 



The estimate by Wood & Reynolds (1999), suggesting that between 5% and 20% of 



the all-sky Ha background may be attributable to scattering by dust appears to be in good 
agreement with our estimates for about half of the high-latitude sky, especially if lower 
latitudes are included (Fig. 4a), but the authors failed to account for the highly structured 
distribution in the ratio of Ha scattering to Ha in-situ emission (Fig. 4b) in their model. It is 
precisely this structure that is important in attempts to determine the small-scale anisotropy 
of the cosmic microwave background at angular scales of < 1 degree. 

Also shown in Fig. 5 are the published positions where WHAM measurements of the [S 



II] /Ha line intensity ratio have been obtained by Reynolds (1985, 1988) and Madsen et al. 
(2006). It appears, the selection criteria applied to these observations have resulted in a 



somewhat biased data set in the sense that the chosen directions fall either very close to the 
Galactic plane or predominantly in high-latitude regions where we identify the scattered Ha 
contribution to be very small. We averaged the ratios of the scattered to total Ha intensities 
predicted by Eq. (6) for all those individual pointings shown as small open circles falling 
outside the masked area in Fig. 5, with a resulting average of of scattered-to-total Ha of 



0.14. Also shown with a cyan outline is a continuous region surveyed by Madsen et al. 



(2006). In Fig. 6 we show the line emission ratios [S II] /Ha as a function of the predicted 



ratio of scattered-to-total Ha for those pointings within this region that fall outside the 
near-plane masked area. For comparison, we also entered the corresponding values for the 
excess emission in LDN 1780, based on our recent WHAM measurements. While the data 



from the Madsen et al. (2006) survey appear to represent to a large extent a typical WIM 



environment with little scattering and significantly enhanced [S II] /Ha ratios compared to 
those seen in H II regions, this comparison shows LDN 1780 to be clearly distinct from the 
WIM, with an emission line ratio similar to that of H II regions. This is an independent 
indication that the Ha excess in LDN 1780 is a result of scattering. Note that the total ratio 
of scattered-to-total Ha emission toward LDN 1780 is about 0.65, reflecting the fact that in 
addition to the Ha excess in LDN 1780 analyzed in this paper, additional Ha intensity of 
both scattering and recombination origin is seen in the direction of LDN1780. 
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It is our conclusion that the ratio of scattered to total Ha background across the high- 
latitude sky is highly variable and that this could be demonstrated by conducting an unbiased 
survey of the [S II] /Ha line intensity ratio. While our results may lead to a reduction in 
the Ha intensities attributable to in-situ emission in the WIM, they do not eliminate in-situ 



emission in the WIM as a real phenomenon requiring further study (e.g. Wood et al. 2010). 



After completion of this work a paper by Seon et al. (arXiv:1006.4419) became avail- 
able, dealing with the analysis of the far-ultraviolet (FUV) continuum background observed 
by SPEAR/FIMS. The latitude-dependent FUV background is consistent with scattering by 
interstellar dust. At Galactic latitudes \b\ > 25° these authors find a close linear correlation 
between the diffuse Ha background taken from the Finkbeiner map and their FUV back- 
ground intensities. The most likely interpretation of this result is that the high-latitude Ha 
background is significantly affected by its dust-scattered component. 



CONCLUSIONS 



We confirm the prediction of Jura (1979) and the findings of Mattila et al. (2007) 



that the excess surface brightness in Ha light of LDN 1780, and by extension in other 
high-latitude dust structures ( Lehtinen et al.pOlO ), is due to scattering of Ha photons 
present in the Galactic ISRF, having originated in emission regions elsewhere. 

Our determination of the line-to-continuum ratio of the scattered light in LDN 1780, 
when combined with earlier work relating the optical continuum intensity to the 100 fim 
surface brightness of optically thin high-latitude clouds, allows us to predict the in- 
tensity of the expected scattered Ha intensity for large fractions of the high-latitude 
Galactic sky. 

We find the ratio of scattered to total observed Ha intensities for many regions domi- 



nated by Galactic cirrus dust to be of the order of 0.5, higher than predicted by Wood 



& Reynolds (1999), but with considerable small-scale structure and spatial variation. 



The standard assumption in the reduction of CMB observations, namely that the ob- 



served Ha all-sky map of Finkbeiner (2003) is a suitable template for free- free emission 



in the Galactic foreground, must be reexamined, especially as it applies to high-latitude 
parts of the sky. 
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Table 1. Filter Characteristics and Exposures 



Filter 


A [nm] 


FWHM [nm] 


Exp. Time [s] 


Obs. Date 


BATC#8 


607 


30.0 


9000 


4 April, 2008 


BATC#9 


666 


51.7 


8100 


4 April, 2008 


BATC#10 


705 


29.0 


9000 


4 April, 2008 


Ha Narrow 


656.3 


6.0 


9000 


11 April, 2008 



Fig. 1. — (top) Red continuum image of LDN 1780, identifying the 378 sampling areas 
(colored circles) and 40 sky reference fields (white). The field shown is 88' x 56'. Each 
sampling area is 74" in diameter. Typical line-of-sight extinctions through LDN 1780 are 
Ay ~ 1 mag (cyan, blue, yellow), Ay ~ 2 mag (green) and Ay ~ 4 mag (red). North is up 
and east is to the left, (middle) 9000 s exposure of LDN 1780 with a 6 nm wide narrow-band 
Ha filter, (bottom) Ha excess image of LDN 1780 after corrections for sky and continuum 
contributions in the narrow-band Ha filter. The image has been smoothed with a 7-pixel 
Gaussian. The intensity scales are counts px _1 . 
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Fig. 2. — Cross correlation of Ha surface brightnesses in LDN 1780, measured on our New 
Mexico Skies (NMS) image vs. the Finkbeiner (2003) map. The symbol colors refer to 
the source areas identified in Fig. l(top). The solid diagonal line indicates the expected 
1:1 correlation; the dashed curve is a linear least-squares fit to the data. On average, the 



Finkbeiner (2003) intensities are about 11% higher. 
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Fig. 3. — (a) Averaged continuum intensities in LDN 1780, measured with the BATC #8 
and BATC #10 filters, are plotted as a function of the optical depth at 160 fim (A v = 1 mag 
corresponds to r(160) /xm) = 0.0006). The symbol colors refer to the different source regions 
identified in Fig. l(top). The solid curve represents the least-squares fit of the function given 
in Eq. (2) to the data, (b) Plot of the excess Ha intensity (NMS data) in LDN 1780 as a 
function of the optical depth at 160 /im . The symbol colors refer to the different source 
regions defined in Fig. 1 (top). The solid curve shows the expected scattered light intensities 
for the case here the Ha line-to-continuum ratio is the same as observed at the sun. The 
dashed curve represents the case where the Ha line intensity is increased by 38% compared 
to the continuum relative to the solar neighborhood. 
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Fig. 4. — (a) Frequency distributions of total observed Ha intensities (black) from the 
Finkbeiner| ( [2003] ) map for the entire high-latitude sky limited by the condition r(Ha) < 1, 
of the predicted scattered Ha component (green) according to Eq. (6), and of the difference 
(red), which we attribute to in-situ recombination emission. — (b) Frequency distributions 
of total observed Ha intensities (black) from the Finkbeiner ( |2003[ ) map for a high-latitude 
region bounded by 1 = 340° and 1 = 60° on longitude, b = +35° and b = +75° in latitude, 
of the predicted scattered Ha component (green) according to Eq. (6), and of the difference 
(red), which we attribute to in-situ recombination emission. — (c) Frequency distributions 
of the total observed Ha intensities before (black) and after (blue) extinction correction 



according to Bennett et al. (2003) compared with the residual Ha intensities after subtraction 



of the scattered light component, again before (red) and after (pink) extinction correction. 
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Fig. 5. — An all-sky map (1° resolution map) of the ratio of the predicted scattered Ha 
intensities to the total observed Ha intensities from the Finkbeiner (2003| map. The blocked- 
out regions of the map are areas where the optically thin condition r(Ha) < 1 is no longer 
valid. The small white circles identify the published positions of WHAM measurements of 
the [S II] /Ha emission line ratio by Reynolds (1985, 1988) and Madsen et al. (2006). The 
region outlined in cyan represents the portion of the sky in which detailed surveys of emission 
line ratios were conducted by Madsen et al. (2006). 
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Fig. 6. — Scatter diagram of individual measurements of the line emission ratio [S II] / Ha 
in the Madsen et al. (2006) survey field identified in Fig. 5 plotted against the predicted 
ratios of scattered-to-total Ha intensity from our l°-resolution map shown in Fig. 5. The 
correspondingly measured values for LDN 1780 are indicated by a star symbol. Note that the 
line emission ratio [S II] / Ha for LDN 1780 is the result of differential spectrophotometry 
characterizing the excess emission in LDN 1780 relative to the surrounding sky regions, while 
the Madsen et al. (2006) survey field data are the result of absolute photometry. 



